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Abstract. The hydrodynamic limit of a one dimensional kinetic 
model describing chemotaxis is investigated. The limit system is 
a conservation law coupled to an elliptic problem for which the 
macroscopic velocity is possibly discontinuous. Therefore, we need 
to work with measure- valued densities. After recalling a blow-up 
result in finite time of regular solutions for the hydrodynamic 
model, we establish a convergence result of the solutions of the ki- 
netic model towards solutions of a problem limit defined thanks to 
the flux. Numerical simulations illustrate this convergence result. 
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1 Introduction 

1.1 Modeling 

Chemotaxis is a process in which a population of cells rearrange its struc- 
tures, reacting to the presence of a chemical substance in the environment. 
In the case of positive chemotaxis, cells migrate towards a concentration gra- 
dient of chemoattractant, allowing them to aggregate. Since several years, 
many attemps for describing chemotaxis from a Partial Differential Equa- 
tions viewpoint have been considered. The population at the macroscopic 
level is described by a coupled system on its density and the chemoattractant 
concentration. The most famous Patlak, Keller and Segel model [181 [23] is 



formed of parabolic or elliptic equations coupled through a drift term. Al- 
though this model has been successfully used to describe aggregation of cells, 
this macroscopic model has several shortcomings, for instance the detailed 
individual movement of cells is not taken into account. 

In the 80's, experimental observations have shown that the motion of 
bacteria (e.g. Escherichia Coli) is due to the alternance of 'runs and tum- 
bles' [T| 114 4 l20 l I22j. Therefore kinetic approaches for chemotaxis have been 
proposed. The so-called Othmer-Dunbar-Alt model [201 1221 [21] describes 
the dynamic of the distribution function / of cells at time t, position x and 
velocity v and of the concentration of chemoattractant S : 

' d t f + v-V x f= [ (T[S](v> v)f{v') - T[S](v v')f{v)) dv' , 

Jv '^ v , (1) 
— AS + S = p(t, x) := / f(t,x,v)dv. 

V Jv&V 

In this equation, V is the set of admissible velocities. The turning kernel 
T[S](v' — > v) denotes the probability of cells to change their velocity from 
v' to v . Several works have been devoted to the mathematical study of this 
kinetic system, under various assumptions on the turning kernel, see for 
instance [11} \10\ \13\ [T7] . Here we shall assume that the velocities of cells 
have the same modulus c, so that V = S c := {v \ \\v\\ = c}. 

Derivation of macroscopic models from (pQ) has been investigated by sev- 
eral authors. When the chemotactic orientation, or taxis, that is the weight 
of the turning kernel, is small compared to the unbiased movement of cells, 
the limit equations are of diffusion or drift-diffusion type. In [161 121]) the 
authors show that the Patlak-Keller-Segel model can be obtained as a dif- 
fusive limit for a given smooth chemoattractant concentration. A rigorous 
proof for the case of a nonlinear coupling to an equation for the chemical 
can be found in [IT], leading to a drift-diffusion equation. 

In this paper we focus on the opposite case, where taxis instead of undi- 
rected movement is dominating. The model has been proposed in [12] , and 
we briefly recall how it is obtained. The limit problem is usually of hyper- 
bolic type, see for instance [15] , Dominant taxis is reflected in the transport 
model by the fact that the dominating part of the turning kernel depends 
on the gradient of the chemoattractant. At this stage, two possible models 
are encountered. On the one hand, we can assume that cells are able to 
compare the present chemical concentrations to previous ones and thus to 
respond to temporal gradients along their paths. The decision to change 
direction and turn or to continue moving depends then on the concentration 
profile of the chemical S along the trajectories of cells. Thus the turning 
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kernel takes the form (independant on v) 

T[S] («'->«) = </>(dtS + v'-V x S). (2) 

On the other hand, if cells are large enough, it can be assumed that they 
are able to sense the gradient of the chemoattractant instantly so that we 
can use instead the expression 

T[S}(v'^v) = 4>(v'-V x S). (3) 

Theoretical results as well as numerical simulations for models (HI)-© and 
CE])-© are proposed in [27]. 

The function (f> in the preceding formulae is the turning rate, obviously it 
has to be positive and monotone. More precisely, for attractive chemotaxis, 
the turning rate is smaller if cells swim in a favorable direction, that is 
dtS + v' ■ V X S > (or v ■ S/ X S > 0). Thus <fi should be a nonincreasing 
function. The converse holds true for repulsive chemotaxis. A simplified 
model for this phenomenon is the following choice for <j>: we fix a positive 
parameter a, and take 

# 60 - w , *'<o, *w = {*; /4 w 

where <fio is a given constant. Moreover, for the sake of clarity of the paper, 
we will assume the following symmetry on <j>: there exists an odd function 
4> such that 

5 ~ ~ ~ 

(j)(x) = <t> (- + <p{x)), and <f)(-x) = -(j>{x). (5) 
o 

The turning kernel ([2]), compared to ([3D, makes drastic changes in the 
behaviour of the solutions to the kinetic model (see |27J). Up to now we 
cannot take it into account in the theory, so that we focus in the follow- 
ing on the expression ([3]). As observed above, this can be considered as a 
biologically relevant model. 

In the turning kernel, a specific parameter quantifies the "memory" of 
the bacteria. When this parameter is small, a specific asymptotic regime 
leads to a macroscopic, hydrodynamic model. In order to introduce this 
parameter, we rescale the system (pQ) by setting 

x = XqX, t = tot, v = VqV, 
S(t,x) = S S(t,x), f{t,x,v) = f J(t,x,v), cf)(z) = <po${z), 
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where 4>q is the typical value for the size of the turning kernel, vq = c is 
the typical speed, xq = L x is the characteristic length of the device and the 
typical time is defined by to = xq/vq. Dropping the bars, the scaled version 
of (H|) reads 

d t f + vV x f = \j {T[S](v' ^ v)f(v>) -T[S](v ^ v')f(v)) dv', 

-AS + S = p, 

where e = ° 1 is the parameter we are interested in: it corresponds 

4> x 

to the time interval of information sampling for the bacteria. The hydro- 
dynamic limit corresponds to e — > 0, and we first recall formally how it is 
obtained. 



1.2 Formal hydrodynamic limit 

We focus in this work on the one dimensional version of ((6l), so that the 
transport takes place in x E R and the set of velocity is V = {— c, c}. The 
expression of the turning kernel simplifies in such a way that ([6]) with (|3|) 
rewrites 

dtfe + vd x f £ = l(<K-vd x S e )f e (-v) - <f>(vd x S e )f e (v)), v G V. (7) 

- d xx S £ + S £ = p £ = f £ (v) + f £ (-v). (8) 
We formally let e go to assuming that S and / admit a Hilbert expansion 

fe = /o + ef x + • • • , S e = S + eSx + ■ ■ ■ 

Multiplying by e and taking e = 0, we find 

(j>(-cd x S )f (-c) = <t>(cd x S )f (c). (9) 

Summing equations (|7|) for c and — c, we obtain : 

dt(fe(c) + U-c)) + cd x (f £ (c) - U-c)) = 0. (10) 

Moreover, from equation we deduce that 

tr\ t( \ tt-cdxSp) -<p{cd x So) 

/o(c) " /o( " c) = H-cd x s o) + Hcd x s ) ifoic) + /o( " c)) 



4 



The density at equilibrium is defined by p := fo(c) + /o(— c) = J fo(v)dv. 
Taking e = in (JTDJ) we finally obtain 



d t p + d x {a(d x S )p) = 0, 

where a is defined by 

f« c ^ (j)(-cd x So) - <j>(cd x S ) 4 
a(9 x 5 ) = c— — - = -c{4>{-cd x S Q ) - <p{cd x S )), 

4>{-cd x S ) + 4>(cd x b ) 5 

and we have used (|5|) for the last identity. Notice that a is actually a 
macroscopic quantity, since we can rewrite 



a(d x S ) = -- [ v<j)(vd x So)dv, 
o JV 



so that this expression is independant of the sign of c. 

We couple this equation with the limit of the elliptic problem ((SJ) for 
the chemoattractant concentration, so that, in summary, and dropping the 
index 0, the formal hydrodynamic limit is the following system 

d tP + d x (a(d x S)p) = 0, (11) 

a(8 x S) = \c (0(-cd x S) - ^{cd x S)), (12) 
5 

-d xx S + S = p, (13) 
complemented with the boundary conditions 

pit = 0, x ) = p mi (x), lim p(t,x) = Q, lim S(t,x) = 0. (14) 

x— >±oo x— >±oo 

The formal hydrodynamic limit from (JT])-© to (|lip -(|12 p -(|13 p has been 
obtained in |12j and proved rigorously in the two-dimensional setting for a 
given smooth S. The aim of this paper is to give an account of the problems 
and open questions arising in the study of the whole coupled system. 



1.3 Preliminary remarks 

First notice that, even in this one dimensional framework, this study 
leads to difficulties mainly due to the lack of uniform estimates for the 
solutions to the kinetic model when e goes to zero and consequently to the 
very weak regularity of the solutions to the limit problem. Even though 
existence of weak solutions to the kinetic model is ensured in a LP setting, 
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no uniform L°° bounds can be expected. The reader is referred to |27j for 
some numerical evidences of this phenomenon, which is the mathematical 
translation of the concentration of bacteria. This is some kind of "blow-up 
in infinite time" , which for e = leads to actual blow-up in finite time, and 
creation of Dirac masses. Moreover the balanced distribution vanishing the 
right hand side of ([7]) depends on S e ; thus the techniques developed e.g. in 
[TT] cannot be applied. 

We turn now to formal considerations about the limit system, noticing 
on the one hand that a solution of (I13p has the explicit expression 



so that the macroscopic conservation equation for p (|lip can be rewritten 



When a is the identity function, this is exactly the so-called aggregation 
equation, which has been studied by several authors, see [21 El SI US] and 
references therein. In particular, finite time blow-up is evidenced when the 
kernel K is not smooth enough. 

On the other hand, taking a = in the definition of P| and assum- 
ing that the chemoattractant concentration is increasing for x < xq and 

decreasing for x > xq (which is usually true when cells aggregate at the 

3 

position xo), we deduce that a(d x S) = — c sgn(x — Xq) which presents a 

5 

singularity at x = xq. The conservation equation (|lip becomes therefore a 
linear conservation equation with a discontinuous compressive velocity field, 
and it is well known that the solution is a Dirac mass. If a is positive, it 
turns out that a Dirac mass appears as well, after a finite time. 

In summary, we have to deal in the limit system with some kind of weakly 
nonlinear conservation equation on the density p. Indeed on the one hand 
the expected velocity field depends on p, but in a nonlocal way. On the 
other hand, this equation behaves like linear equations with discontinuous 
coefficients, in the sense that it admits measure- valued solutions. Therefore 
a major difficulty in this study will be to define properly the velocity field 
a = a(d x S) and the product ap. 

The paper is organized as follows. In Section[2]we consider the aggregation- 
like equation (116p . and recall existence and uniqueness results as well as the 
existence of a finite time for which L°°-weak solutions of (jlip -(|13 p blow up. 
In Section [3l we investigate the hydrodynamical limit of system ©-([8]) and 
prove in particular that it gives rise to a somehow natural definition of the 



S(t, x) = K* p(t, .)(x), where K(x) = \e~\ x \ 



(15) 



dtp + d x (a(d x K * p)p) = 0. 



(16) 
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flux in the conservation equation. Some numerical simulations illustrating 
this result are furnished in Section [U Finally, we end this work with some 
conclusions and remarks. 



2 Aggregation-like equation 

In this section, we consider the equation 

/ d t p + d x {a{d x K * p)p) = 0, , . 

1 p(t = 0,-)=p mi . [U) 

where K is given by (fT5|) . We assume that 

< p mi £ L 1 nL°°(M). (18) 

When a(x) = x, this equation is the so-called aggregation equation (see e.g. 
[21 HI IH]). It is known that for singular d x K, solutions blow up in 
finite time. More precisely, we show the blow-up in finite time of L°° weak 
solutions. Most of the results presented in this section are obtained thanks 
to a straightforward adaptation of techniques developed in [21 O [6], [19]. 
Therefore some proofs are not detailed. 



2.1 Existence and uniqueness of local L^-weak solution 

We prove in this section the local existence and uniqueness of a solution. 

Theorem 2.1. Let p ini € L 1 nL°°(R). Then there exists a T > such that 
there exists a unique weak solution p to fj 1 7f) / moreover p E C([0, T]; L 1 n 
L°°(R)). 

The proof is an adaptation of results in [21 [31 [19]. We first recall the 
definition of the characteristics for this system: X(s; x, t) is a solution of the 
ODE 

— (s;x,t) = a(d x K*p){s,X{s)), X(t;x,t) = x. (19) 
as 

Then we have the following representation of the solution of (j 1 T p : 

p(t,x) = p im (X{0;x,t))exp (- J d x a(d x K * p)(s, X(s - t;x,t)) ds) . 

(20) 

The proof of this theorem relies strongly on the following estimates: 
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Proposition 2.1. Let p ini such as in (|18[) and let p be a solution of ()17|) 
on [0, T]. Then there exists T > suc/i f/mf for all t G [0, T], f/iere exists a 
nonnegative constant C such that 

\\p(t, OIIlKR) + IM*> ')IU°°(K) < c, 
where C only depends on Wp^Wl 1 ^) an d ||p im |U°°(K)- 

Proof. The L 1 estimate is an easy consequence of the mass conservation. 
Then, 

\d t p + a(d x K * p)d x p\ = \ - d x (a(d x K * p))p\ < 2||a'|| 0O |p| 2 . 
Integrating along the characteristics curves, we get 

||p(V)Hl» < ||p m i|L- + 2||a'|| 00 f \\p(s, floods. 

Jo 

We deduce that as long as 2||a'||oo||p im ||L°°i < 1, 

||P(V)IU°° < --, ^TTTii [T^JTl 7- 

We notice that T should satisfies the bound T < l/(2||a / || 00 ||p m2 ||^oo). □ 

Proof of existence. We do not detail the proof of the existence of solution 
which can be deduced thanks to an adaptation of [21 El 02], where the study 
of an aggregation equation is proposed. We just recall the main argument 
of the proof in the following steps : 

1. We construct a family of approximating solutions (p e ) by solving (fTT|) 
with initial data p mi * g £ where g e is a mollifier. 

2. We state uniform Lipschitz estimates in space and time on the se- 
quences (a(d x K * p e ))e and (X e ) e and use the Arzela-Ascoli Theorem 
to extract converging subsequence. 

3. We pass to the limit in the representation (|20l) . □ 

Proof of uniqueness. The idea of this proof is to use the quantity S. 
Since this idea will be developed for measure- valued solutions, we detail 
this proof. Computations are done for regular solutions, nevertheless they 
can be made rigorous by introducing a regularization and passing to the 



S 



limit (see [2]). Let us consider two classical solutions p\ and p 2 . Denoting 
ai = a(d x K * pi) for i = 1, 2, we have 

dt(pi ~ P2) + d x (ai(p! - p 2 )) + d x {p 2 {ai - a 2 )) = 0. (21) 

Define S(t,x) := (d x K * (pi — p 2 )(t, -))( x ) which solves the problem 

-d xx S + S = Pl -p 2 , onl. (22) 

We notice that when t = 0, we have <S(0, x) = 0. From the weak formulation 
of equation ([2T]) with the test function S 1 , we have 



where 



f [ dt(pi - p 2 )S dxds = I + II 
Jo Jr 

1=1 / ai(pi - p 2 )S dxds, 
Jo Jr 

11= / p 2 (ai — a 2 )S dxds. 
Jo Jr 



(23) 



For the term I, we have using (J22J) and integration by parts 

I=[ [ ai(—d xx S+S) S dxds = - I I d x ai\d x S\ 2 dxds+ I f a\S 2 dxds, 
Jo Jr 2 Jo 7r jo Jr 

Moreover, 

4 8 
d x ai = --c{<t> {-cd x S\) + (j) (cd x Si))d xx Si < max{-c||0 || £oo 5i,0}, 

where we use the fact that is a nonincreasing positive function. From 
the L°°-bound on Si, we deduce that there exists j3 € L 1 ([0,T]) such that 
d x ai < (3. Thus 

I<\l [ \d x S\ 2 dxds+ [ [ ai S 2 dxds. 
2 jo Jr Jo Jr 

Then, the estimate |ai| < |c gives 

I< \ [ l3(s)\\d x S\\ 2 L2 ds + lc f \\S\\ 2 L2 ds. (24) 
1 Jo " Jo 
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For the term II of ()23|) . we have thanks to the Cauchy-Schwarz inequality 



I -Til 



/ / p 2 (a(d x K * pi) - a(d x K * p 2 ))S dxds 
Jo Jr 

/ ||/?2 ||l«"||o'||do||^x^ * (Pi - P2)\\lA\ S \\l 2 ds 

Jo 



< 

Since d x S = d x K * (pi — p 2 ), we obtain 



i r* 

\H\ < 77 / llPzlU-llo'lloodl^^Hia + \\S\\l 2 )ds. 

1 Jo 



(25) 



Then, we notice that using (|22p and thanks to an integration by parts the 
left hand side of (1251) can be rewritten 



f ! dt(pi - p 2 )S dxds = \ [ (\d x S\ 2 + S 2 )dx (26) 
Jo Jr 1 Jr 

Finally, we deduce from (p|) . (plj) . (l25l) and (|26l) . 

^(|^5| 2 + S 2 )dx < J ||p2|U«||a1oo (p(s)\\d x S\\ 2 L2 + ~c\\S\\ 2 L ^ ds. 
Uniqueness follows from a Gronwall type argument. □ 

2.2 Blow-up in finite time 

The blow-up of solutions of a one dimensional aggregation solution is pro- 
posed for instance in [6] where it is proved by the method of characteristics 
that aggregation of mass occurs. In [2_1[3], the finite time blow-up is obtained 
thanks to an energy estimate. We assume that the initial data is given sym- 
metric with respect to and positive. It is easy to show then that for all 
t > 0, p(t,x) = p(t,—x). Moreover, for the sake of simplicity, we assume 
that there exists 5 > such that supp(p mi ) C [—5, 5]. Then for all x > 5 the 
function S im = K * p im satisfies d x S mi (x) < 0, so that the characteristics 
defined by (fT9|) are inward. Thus for all t > 0, supp(p(t, •)) C [—5, 5]. 

The energy of the system is defined as 

E(t) = I ! {\d x S\ 2 + \S\ 2 ) dx = \ [ pSdx, (27) 
^ Jr * Jr 

where the last formulation is obtained by integration by parts. On the one 
hand, we have the obvious bound 

E{t) < ±\\p\\ L i\\K * p\\ Loa < ^\\pf Ll . (28) 
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On the other hand, using (|17p we have 



^-E(t) = [ a(d x S)d x Spdx. 
at 

Moreover ||<9 a; »S'||_L°o < |||/o tn *||i,i- Since the function a is assumed to be regu- 
lar, there exists ( > such that a(x)x > (\x\ 2 for all x € [— |||p tm ||ii, lllp^lli 1 ]- 
Thus, 

4e(*)>C / \d x S\ 2 pdx. (29) 
at J R 

We now make use of the following result whose proof is given in [3] : 

Proposition 2.2. There exists a constant C > such that for all 5 suffi- 
ciently small, we have for any symmetric nonnegative function p in L 1 (M) 
with a compact support in [—5,5], 

[ \d x K*p\ 2 pdx > C. 
Jm. 

Then, from (|29|) there exists a constant C > such that for all t > 0, 

E(t) - E(0) > Ct. 
Therefore, with (j28|) we have proved 

Theorem 2.2. Let p be a symmetric solution of ()17|) with symmetric, 
positive initial data with compact support included in [—5, 5]. For sufficiently 
small 5, there exists a time T* > for which the solution p ceases to exist, 
i.e. 

I™ l|p(*,-)lliP(K) = forpe (l,oo). 



3 Convergence for the kinetic model 

In this section we investigate the convergence of a sequence of solutions to 
the microscopic model ([7|)-(j8]). We are not able yet to obtain rigorously (fTTI) - 
(fT3j) . We actually prove that the whole sequence of solutions is convergent, 
and that the macroscopic density satisfies a conservation equation with a 
uniquely determined flux. More precisely, the main result of this section is 
the following theorem. We introduce the macroscopic densities 

Pe= fedv, p= fdv. 
Jv Jv 
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Theorem 3.1. Let T > and let us assume that p mi is given in Mbffi)- 
Let (f e ,S £ ) be a solution to the kinetic-elliptic equation ©-([8]) with initial 
data fl m such that p l ™ := f v /* m dv = n e * p m% where r\ e is a mollifier. Then 
as e — > 0, the sequence converges to (p, S) in the following sense : 

p e ^p in S M :=C([0,T];M b (R)-a(M b ,C Q )), 
S £ ^S in CQO,^;!^ 1 ' 00 ^)) -weak, 

and (p, S) is the unique solution in the distribution sense of 

J d tP + d x j = o, 

\ -d xx S + S = p, 
complemented with initial data p mi and where 

J = -d x (A(d x S)) + a(d x S)S a.e. 



(30) 



Before turning to the proof of this result we notice that the problem ()30|) 
is equivalent to 

8 t S - d x K * [d x {A{d x S)) + a(d x S)S] = 0, in V'(R). (31) 

This is obtained by taking the convolution with K of the first equation in 
(|30p . This emphasizes the key role of S in the study of the limit. 



3.1 Preliminary results 

First we recall the following statement on the kinetic-elliptic problem. 

Theorem 3.2. Let T > and e > 0. Assume fi ni € C(M). Then problem 
complemented with initial data /* m admits a unique weak solution 
in C([0,T] x K x V) x C([0, T];C 2 (R)). Moreover, we have the following 
estimates uniform in e > : 

\v\ k f £ dxdv = \v\ k \p ini \(R), ken. (32) 

v 

Proof. The proof of the existence can be found in [27] . The estimates (1321) 
rely on the conservation of the mass and on the fact that since v E V = S c , 
|f| is constant. □ 

Then, we furnish a convergence result for a sequence of functions S. 
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Lemma 3.1. Let (/O ra ) n( =N be a sequence of measures that converges weakly 
towards p in Sm- Let S n (t,x) = (K * p n (t,-))(x) and S(t,x) = (K * 
p(t, -))(x), where K is defined in (I15p . Then when n — > oo we have 

d x S n {t,x) — > d x S(t,x) for a.e. t G [0,T], i£l, 
d x S n (t,x) — 1 d x S(t,x) in L°°w — * 



Proof. We have that 

d x S n (t, x) = (d x K * p n (t, •)) (x) 



1 x — y 
2\x — y\ 



-\ x -y\p n (t,dy). (33) 



Let £ > 0, we regularize the convolution kernel by introducing the following 
functions : 



<t>x,s{y) 



\t~ — H" e '* y '> on (— oo,x — e] U [x, +oo), 
2 \x — y\ 

\\lh {y ~ x) + l )' on(x-e,x). 



1 x — y 



\ x ~y\ 



2 \x — y\ 
1 fl + e~ £ 



on (— oo, x] U [x + e, +oo) 



(y — x) — 1 ) , on (x, x + e). 



2 \ £ 

With this definition, we clearly have for all x, y in 

1 x — y 



2\x — y 

Moreover by definition of the weak convergence, 



lx - yl < <t> x Jv). 



(34) 



lim / 4> X £ (y) p n (t,dy) = / <f> X £ (y)p(t,dy). 

n— >+oo 



Then, from (|33l) and flSj), we deduce 



limsupd x S n < / <j>x,e{y) p(t,dy). 

n->+oo JR 



Moreover, 

/ <f>x,e(y) p(^ d y) 

Jr 



d x S + 



1 



1 



1 + e~ e 
1 H (y - x) 



<^5+-(l-e- £ )|p(t, 
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By the same token with ift x ,e, we obtain the estimate 

d x S e- £ )\p(t, -)|(R) < Hminf d x S n 

< limsupdA < 8 X S + 1(1 - e- e )\p(t, .)|(R). 

n— >+oo * 

Letting e — > 0, we get lim n _j. +00 d x S n (t, x) = d x S(t,x) for almost all t € 
[0, T] and iGl. □ 
We turn now to the uniqueness of S, which is the key point to get the 
uniqueness result in Theorem 13.11 



3.2 Uniqueness for 5* 

In [2] (see proof of Theorem 12.11 too) . the authors obtain the uniqueness 
on the aggregation equation by introducing a quantity which appears here 
naturally to be the potential S. They get an estimate that relies strongly on 
the L°° bound on the density whereas here it is only measure- valued with a 
finite total variation. Therefore, we have to work in a weaker space, and we 
use the fact that the function S defined by S = K * p is a weak solution of 
(I3ip . We have the following result. 

Proposition 3.1. Let Si and S 2 be two weak solutions of ([M]) in L°°([0, T]; B(R))n 
C([0,T]; ^^(I)) with initial data S{ ni and Sf™ respectively. Then there 
exists a nonnegative constant C such that 

|Pl — '-'2||l° o ([0,T];W 1 . — ^ H^l _ *2 II W 1 - 1 ^) ■ 

Proof. We introduce a function A E C°°(R) such that A 1 = a. Now, 
differentiating (j3T|) and noticing that K satisfies —d xx K + K = 6q, we get 

d t d x S + d x {A(d x S)) - d x K * A(d x S) + K * (a(d x S)S) - a(d x S)S = 0. (35) 

The definition of S, S(t,x) = (K * p(t, -))(x), implies that d x S belongs to 
L°°(0,T; BV(M)). Therefore equations f)31j) — (|35|) have a sense in their weak 
formulation. Let S± and 52 satisfy the weak formulations of (|31|) — f)35 |) with 
initial data S\ m and S™ 1 respectively. We denote by a% = a(d x S\) and 
a2 = a(d x S2)- We deduce from ([35]) that 

dtd^S! - s 2 ) + d x (A(d x Si) - A(d x s 2 )) = 

d x K * (A(d x Si) - A{d x S 2 )) + a x S x - a 2 S 2 - K * ( ai Si - a 2 S 2 ). 
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Multiplying this equation by sign (c^ (Si — S2)), integrating with respect to 
x and using the properties of the convolution product, we deduce 



d_ 

dt 



[ \d x (Si - S 2 )\dx < II^KHoo f \A(d x S 1 )-A(d x S 2 )\dx+ 

JR JR 

+(1 + Hi^lloo) / |oi5i — a 2 52| dx. 

Jr 



The function a being regular, we deduce 
d 



dt 



[ \d x (Si-S 2 )\dx<C [ \d x (S l -S 2 )\dx + C l [ \S 1 -S 2 \dx. (36) 
JR Jr Jr 

By the same token with equation f)31 1) . it leads to 

f \Si - S 2 \ dx<C 2 ( \d x (Si - S 2 )\dx + C 3 ( \S X - S 2 \dx. (37) 
Jr Jr Jr 



d_ 

dt 



Summing ([37]) and ([36]) . we deduce that there exists a nonnegative constant 
C such that 

— US'! - 52||w 1 - 1 ()R) < C\\Si - S 2 \\ w i,i(^y 

Applying the Gronwall Lemma allows to conclude the proof. □ 



3.3 Proof of Theorem EH 

Let (f £ ,S £ ) be a solution of ©-([8]). For fixed e > 0, we have f £ € 

C([0, T] x IR x V). We define the flux J £ := / vf £ dv and the macroscopic 

Jv 

velocity 

in q \ j v v(t>{vd x S £ )dv 4 
a{d x S £ ) = — f = -v{(j){-vd x S £ ) - 4>{vd x S £ )). 

Jv 9\vd x b £ ) dv 5 

We can rewrite the kinetic equation (|7|) as 

1 5 

dtfe + vd x f £ = -(<p{-vd x S £ )p £ - -f £ ). 

Taking the zeroth and first order moments, we get 

d tPe + d x j £ = 0, (38) 

8 t J £ + v 2 8 x p £ = ^( a (d x S £ )p £ - J £ ). (39) 
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From ([38]) . we deduce that 



Vte[0,T], \p £ (t,-)\(R) = \p mi 

Therefore, for all t S [0, T] the sequence (p e (t,-)) e is relatively compact in 
M b (R) - a{M b (R),C {R)). Moreover, there exists u £ e L°°([0, T],BV(R)) 
such that p e = d x u £ . From ([38|h we get that dtu £ = — J £ and with estimate 
(j32|) for A; = 1 we deduce that u £ is bounded in Lip([0, T], It implies 

the equicontinuity in t of (p £ ) £ . Thus the sequence (p £ ) £ is relatively compact 
in Sm and we can extract a subsequence still denoted (/5 e ) £ that converges 
towards p in Sm- 

We recall that S e (i, x) = (K*p £ (t, -))(x) where K(x) = ^e~^. Denoting 
S(t,x) := (K * p(t,-))(x), since p G S M , we have 5 G L°°([0, T]; J3V(M)). 
From Lemma l3.1| the sequence (d x S £ ) £ converges in L°°w — * and a.e. to 
d x S as e goes to 0. 

From (|38jH(|39p . we have in the distribution sense 



dtPe + d x (a(d x S £ )p £ ) = d x {a(d x S £ )p £ - J £ ) = ^ed x (d t J £ + v 2 d x p £ ) = R £ . 

(40) 

Now, for all V £ C c 2 ((0,T) x R), we deduce from ([32D 

{d t J £ + v 2 d x p £ )d x ^dxdt < |«|| /0 i ™|(lR)||^9 :c V||L- + ^| 2 |p im |(lR)||5 :E:E Vl|L- 

This implies that the limit in the distribution sense of the right-hand side 
R £ of (j4"0|) vanishes. 

On the one hand, multiplying equation (|8|) by a(d x S £ ) and introducing 
the real-valued function A such that A' = a, we get 

a{d x S £ )p £ = -d x (A(d x S £ )) + a(d x S £ )S £ , (41) 

so that we can rewrite the conservation equation (|40p as follows, in T>'(R) : 

dtPe + d x {-d x A{d x S £ ) + a(d x S £ )S £ ) = \ed x (d t J £ + v 2 d x p £ ). (42) 

5 



Taking the limit e — > in the distribution sense of equation (J42J) , we get 
that in V'(R) 

dtp + d x {-d x A{d x S) + a(d x S)S) = 0, (43) 

where S(t,x) = {K * p(t, -))(x). We recall that we have chosen the initial 
data such that p l £ % = rj £ * p mi where r\ £ is a mollifier. Therefore p % £ % — v p m ' 
mM b {R)-a{M b (R),C (R)). 
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On the other hand, as noticed above, S satisfies (|3ip and f)35[) in the 
distribution sense. Proposition 13. II above asserts that S satisfying (f3Tj) - ([35"j) 
is unique. Thus p is unique since, if we assume that there exist p\ and 
P2 satisfying ([4~3]) in the distribution sense, then by the uniqueness of the 
solution of (|31|) ~(|35 p . we have that K * p\ = K * p2 which implies that 
p\ = P2- Finally, thanks to the uniqueness, all the sequence p £ converges to 
P in S M - 

4 Numerical simulations 

We illustrate the previous convergence result with some numerical sim- 
ulations of the problem ©-([8]). We discretize the kinetic equation thanks 
to a semi-lagrangian scheme and the elliptic equation for S is discretized 
with P\ finite elements. We refer the reader to |27j for more details on the 
numerical scheme. Notice that letting e go to in the simulations is very 
difficult because of the high numerical diffusivity of the scheme. 

We have chosen to present simulations with realistic numerical values. 
For the bacteria Escherichia Coli the velocity is c = 20. 10 -6 m.s~ l and the 
density of cells is no = 10 11 m _1 . The domain is assumed to be an interval 
of length xo = lcm. The turning kernel is given by (j3|) with (f> in (j4]). Due 
to the large value of no, the value of the parameter a should be very large to 
have an influence; thus this parameter does not play a role in the dynamics 
of bacteria and for the simulations we have fixed a = 1. We assume that 
the initial concentration of cells is a Gaussian centered in the middle of the 
domain. We run simulations with three different values for fa : 4>o = 0.05, 
1 and 20 so that e = vq/^qXq) takes the values 10~ 4 , 2. 10~ 3 and 4. 10~ 2 . 

In Figures [Hand [2] we present evolution of the density of cells with respect 
to the time and to e. We observe the aggregation of cells in the center of 
the domain which is the first step of the formation of a Dirac. As e — > 0, 
the aggregation phenomenon is faster and the solution seems to converge to 
a Dirac. We display the evolution of the gradient of the chemoattractant 
concentration d x S in Figures [3] and [U A singularity in the center of the 
domain appears clearly. 

5 Conclusion 

In this work we have studied the convergence of a kinetic model of cells 
aggregation by chemotaxis towards a hydrodynamic model which appears 
to be a conservation law coupled to an elliptic equation. Although the limit 
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0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009 

x(m) 



0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.006 0.009 0.01 

x (m) 



Figure 1: Time evolution of the density p of bacteria for different values of 
the parameters e. Left : e = 4. 10~ 2 . Right : e = 10 -4 . 




0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009 0.01 



x(m) 

Figure 2: Density p of cells for different values of the parameters e at time 
t = 80 s. As e becomes smaller the concentration effect is more important. 



of the macroscopic quantity p £ and S e have been obtained in Theorem I3.ll 
this mathematical result is not completely satisfactory since the limit model 
(|30p does not allow to define a macroscopic velocity for the flux. Formally, 
this macroscopic velocity is given by a(d x S) defined by (fT2|) . However, since 
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-2 ' ' ' 1 ' 1 1 ' 1 1 

0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009 0.01 

x (m) 

Figure 3: Time dynamics of the gradient of potential d x S with e = 10~ 4 
As time increases the derivative of the potential tends to become singular. 



p is only measure-valued, d x S belongs to BV(U), hence we cannot give a 
sense to the product a(d x S)p. 

A possible convenient setting to overcome this difficulty is the notion of 
duality solutions, introduced by Bouchut and James [7]. In this framework, 
we can solve the Cauchy problem for conservation equations in one dimen- 
sion with a coefficient a that satisfies a one-sided Lipschitz condition. The 
theory in higher dimensions is not complete [9], and Poupaud and Rascle 
|26j (see also[5]) for another approach, which coincides with duality in the 
1-d case. It is actually not difficult to prove that a defined in (fT2|) is one- 
sided Lipschitz. In fact, from p > 0, we deduce that —d xx S < S. After 
straightforward computation, we get 

d x (a(d x S)) = ~c((/>'(-cd x S) + </>'(cd x S))d xx S. 
5 

Therefore, 4> being nonincreasing and smooth, we deduce 

d x {a{d x S)) < max{|c||0 / || Looi S,O}. 
5 
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0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009 0.01 
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Figure 4: Gradient of the potential d x S for different values for e at time 
t = 80 s. 



And the properties of the convolution lead to 

\\S(t,-)\\ L ^ <^\p(t,-)\(R) = ^\p mi \(R). 

Finally, a satisfies the OSL condition : 

3/3 G L 1 ([0, T]), d x a(t, •) < (3(t) in the distribution sense. 

However, we are not able to prove the uniqueness of the duality solutions 
for the hydrodynamic problem. In fact, the uniqueness proof in Section [3. 21 
relies on the fact that the potential S satisfies equation (I3ip and thus on 
the definition the flux J in (|30p . In the framework of duality solution, the 
conservation equation is not a priori satisfied in the distribution sense. A 
generalized flux that has a priori no link with J in (|30p is then introduced. 
The relation between these flux and therefore the passage from J to the 
macroscopic velocity a is still an open question. 

Acknowledgments. The authors acknowledge Benoit Perthame for driv- 
ing their attention on this problem of hydrodynamic limit. 
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